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We report quantum Monte Carlo (QMC), plane- wave density- functional theory (DFT), and inter- 
atomic pair-potential calculations of the zero-temperature equation of state (EOS) of solid neon. We 
find that the DFT EOS depends strongly on the choice of exchange-correlation functional, whereas 
the QMC EOS is extremely close to both the experimental EOS and the EOS obtained using the best 
semiempirical pair potential in the literature. This suggests that QMC is able to give an accurate 
\ treatment of van der Waals forces in real materials, unlike DFT. We calculate the QMC EOS up 

to very high densities, beyond the range of values for which experimental data are currently avail- 
able. At high densities the QMC EOS is more accurate than the pair-potential EOS. We generate 
a different pair potential for neon by a direct evaluation of the QMC energy as a function of the 
separation of an isolated pair of neon atoms. The resulting pair potential reproduces the EOS more 
accurately than the equivalent potential generated using the coupled-cluster CCSD(T) method. 



oo 

o 



c/5 



OO 

o 



- 1—^ 

X 
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INTRODUCTION 



One of the most important goals of ab initio computational electronic-structure theory is the development of 
accurate methods for describing interatomic bonding. Quantum Monte Carlo (QMC) techniques are useful in this 
regard as they can provide a highly accurate description of electron correlation effects. Although QMC methods are 
computationally expensive, they can be applied to systems which are large enough to model condensed matter. 
, In this study we have considered solid neon, in which the bonding arises from the competition between short-range 
repulsion and the van der Waals attraction between the atoms. High-quality experimental measurements of the 
equation of state (EOS) exist, which can be used as reference data. Furthermore, various neon pair potentials have 
O . been developed using experimental and theoretical data, which can also be used for comparison purposes. Solid neon 
is therefore an ideal system in which to test the descriptions of van der Waals bonding and short-range repulsion 
7—1 . offered by various theoretical methods. 

' We have calculated theoretical EOS's for crystalline neon using the QMC and density-functional theory (DFT) ab 
initio electronic-structure methods as well as various interatomic pair potentials. Standard DFT methods do not 
describe van der Waals bonding accurately, but they might be expected to work quite well at high densities, where 
the short-range repulsion dominates. The high-pressure properties of neon are of some experimental interest, because 
neon is often used as a pressure-conducting medium in diamond-anvil-cell experiments i* We have therefore extended 
the range of our QMC EOS for neon to very high pressures (about 400 GPa). 

The zero-point energy (ZPE) of the lattice- vibration modes makes a small but important contribution to the total 
energy of solid neon. We have therefore studied the lattice dynamics of solid neon within the quasiharmonic-phonon 
approximation and within the Einstein approximation, using DFT methods and pair potentials. 

For some time there has been considerable interest in developing neon pair potentials in order to test the accuracy 
of theoretical methods for calculating the properties of materials^ We have performed a direct calculation of the neon 
pair potential using QMC. We compare the accuracy of the EOS predicted by this pair potential with the results 
obtained using other pair potentials, including one obtained from coupled-cluster CCSD(T) calculations^ 

Detailed information about our computational methodologies is given in Sec. [TTJ and DFT calculations of the phase 
stability and band gap of solid neon are reported in Sec. IIII1 The calculation of a neon pair potential using QMC is 
described in Sec. IIVI The lattice dynamics of solid neon are studied in Sec.fV] We compare the EOS's obtained using 
different methods in Sec. IVIl Finally, we draw our conclusions in Sec. IVIII 

Hartree atomic units (a.u.) are used throughout, in which the Dirac constant, the magnitude of the electronic 
charge, the electronic mass, and 4-7T times the permittivity of free space are unity: H = |e| = m e = Aire^ = 1. 
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II. METHODOLOGY 

A. DFT calculations 

1. DFT total-energy calculations 

Our DFT calculations were performed using the CASTEP plane-wave-basis code^ The local-density approxima- 
tion (LDA) and Perdew-Burke-Ernzerhof (PBE) generalized-gradient-approximation 5 exchange-correlation function- 
als were used. The Ne 8+ ionic cores were represented by ultrasoft pseudopotentials* 4 - The EOS calculations were 
performed using a 4 x 4 x 4 Monkhorst-Pack k-point mesh and a plane- wave cutoff energy of 200 a.u., for which the 
DFT energies have converged to about 7 significant figures. The self-consistent- field calculations were judged to have 
converged when the fractional change in the energy was less than 10~ n . The DFT band-gap calculations reported in 
Sec. IIIIBl were performed using the same parameters, except that the plane-wave cutoff energies ranged from 100 a.u. 
for the lowest densities to 800 a.u. for the highest densities. 

2. DFT force- constant calculations 

We used the quasiharmonic approximation^ to evaluate the DFT ZPE of the lattice- vibration modes and we used the 
method of finite displacements and the Hcllmann-Feynman theorem to evaluate the density-dependent force constants. 
Symmetry and Newton's third law were imposed iteratively on the matrix of force constants.™ The DFT force-constant 
calculations were carried out using a plane- wave cutoff energy of 60 a.u., a 3 x 3 x 3 Monkhorst-Pack k-point mesh, 
and ultrasoft pseudopotentials^ The force constants were converged to about 10~ 6 a.u. with respect to the plane- wave 
cutoff energy and the k-point mesh. In the production force-constant calculations, the displacement of the neon atom 
from its equilibrium position was 2.12% of the nearest-neighbor distance in each case, which ensures that anharmonic 
effects are negligible. The force-constant calculations were carried out in both 2x2x2 and 3x3x3 supercells of the 
primitive unit cell, and the difference in the resulting ZPE's was found to be negligible. The dispersion curves shown 
in Sec. W\ were produced using a 3 x 3 x 3 supercell, while the ZPE's that were combined with the static-lattice EOS's 
were calculated in a 2 x 2 x 2 supercell. 

3. DFT orbital-generation calculations 

DFT-LDA calculations were performed in order to generate orbitals for the trial wave functions used in the QMC 
calculations. The QMC calculations made use of relativistic Hartree-Fock neon pseudopotentials^! 9 - and these were 
also used in the DFT orbital-generation calculations. The Hartree-Fock pseudopotentials are much harder than the 
ultrasoft pseudopotentials. Plane-wave cutoffs in excess of 250 a.u. were used in each orbital-generation calculation, 
so the DFT energy was converged to around 10 -3 a.u. This cutoff is very large by the normal standards of DFT 
calculations, but there is evidence that using large basis sets reduces the variance of the energy in QMC calculations. 10 

B. QMC calculations 

1. VMC and DMC methods 

In the variational quantum Monte Carlo (VMC) method, expectation values are calculated using an approximate 
trial wave function, the integrals being performed by a Monte Carlo technique. In diffusion quantum Monte Carlo 11 ! 12 
(DMC) the imaginary-time Schrodinger equation is used to evolve an ensemble of electronic configurations towards 
the ground state. The fermionic symmetry is maintained by the fixed-node approximation,— in which the nodal 
surface of the wave function is constrained to equal that of a trial wave function. Furthermore, the use of nonlocal 
pseudopotentials to represent the Ne 8+ cores necessitates the use of the locality approximation, 14 which leads to errors 
that are second order in the quality of the trial wave function^ 

Our QMC calculations were performed using the CASINO code^ The trial wave functions were of Slater-Jastrow 
form, with the orbitals in the Slater wave function being taken from DFT calculations and the free parameters in 
the Jastrow factor being optimized by minimizing the unreweighted variance of the energy. 17 ' 18 The DFT-generated 
orbitals were represented numerically using splines on a grid in real space rather than an expansion in plane waves in 
order to improve the scaling of the QMC calculations with system size . 19 ' 20 The Jastrow factors consisted of isotropic 



3 



electron-electron, electron-nucleus, and electron-electron-nucleus termsi^l The electron-electron terms describe long- 
ranged correlations and therefore play the most important role in describing van der Waals forces. 



2. Finite-size bias 



The QMC simulations of crystalline neon were carried out in supercells of finite size subject to periodic boundary 
conditions. The electrostatic energy of each electron configuration was calculated using the Ewald method. 22 The 
QMC energy per atom obtained in a finite cell differs from the energy per atom of the infinite crystal due to single- 
particle finite-size effects and Coulomb finite-size effects. The former result from the fact that the allowed k points for 
the Bloch orbitals form a discrete lattice, so that the single-particle energy components change when the size of the 
simulation supercell is changed. The latter, which are the more important in insulators, are caused by the interaction 
of the charged particles with their periodic images. At any given instant, each electron feels itself to be part of an 
infinite crystal of electrons ! 23 ! 24 The resulting bias is negative, and is generally believed to fall off as 1/N, where N 
is the number of atoms in the simulation cellj 11 ! 25 ' 26 ! 27 

In order to eliminate the finite-size bias, simulations were carried out in supercells consisting of 3 x 3 x 3 and 4x4x4 
primitive unit cells. The error in the DFT results arising from the use of a 3 x 3 x 3 k-point mesh is small (about 
0.0001 a.u.), so we conclude that single-particle finite-size effects are negligible. The assumed form of the Coulomb 
finite-size bias was therefore used to extrapolate the results to infinite system size. The static-lattice energy per atom 
in the infinite-system limit is given by 

E s J:(V) = E^(V)+ b -^, (1) 

where Ef^(V) is the Vinet fit (see Sec. Ill D| ) to the DMC static- lattice energy- volume data obtained in a set of iV-atom 
simulation supercells, V is the primitive-cell volume, and b(V) is a parameter determined by fitting. Since we only 
have energy-volume data for two different system sizes, N and M, we may eliminate b(V) and write 

SL NE^{V)-MEl}(V) 
E ™ {V) = N~M ■ (2) 

The pressure due to the static-lattice energy at infinite system size is given by 

J3L rin _ dE% _ Npf{V)-Mp^(V) 

where Pj^(V) = —dEf^/dV is the static-lattice pressure in an iV-atom simulation supercell. 

The zero-temperature, static-lattice energy-volume curves of neon, calculated using DMC in different sizes of sim- 
ulation supercell are shown in Figs. [T] and [2] Vinet EOS's are fitted to the data. The corresponding pressure- volume 
data are shown in Figs. [3] and |4] It can be seen that the DMC pressure- volume curves converge steadily with system 
size, and that the DMC pressure extrapolated to infinite system size using Eq. © is close to the pressure of the 
4x4x4 supercell. This implies that the error introduced by the extrapolation is small, because the extrapolation is 
itself a small correction. 



3. Time-step bias 



The fixed-node DMC Green's function is only exact in the limit of zero time step; the use of a nonzero time step 
biases the DMC energy. An example of the bias in the DMC energy of a pseudoneon crystal is shown in Fig. [5] 
On the other hand, as shown by the 2x2x2 supercell results in Figs. [3] anddj the pressure is very insensitive to 
the time step. We used time steps of 0.005 a.u. and 0.02 a.u. in our production calculations for the 3x3x3 and 
4x4x4 supercells, respectively. A time-step of 0.002 a.u. was used in the DMC pair-potential calculations. The 
target population was at least 320 configurations in each case, while a target population of 1000 configurations was 
used for the pair-potential-generation calculations. 
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FIG. 1: (Color online) Low-density static-lattice DMC energy of FCC neon as a function of volume, evaluated in simulation 
supercells consisting ofnxnxn primitive unit cells using different time steps r. 



C. Pair-potential calculations 

1. Forms of neon pair potential 

We have used the following forms of pair potential: (i) the HFD-B potential proposed by Aziz and Chen^S with 
the parameter values given by Aziz and Slama n 29 i 30 ; (ii) the form of potential proposed for helium by Korona et 
al.r^ containing the parameter values determined by Cybulski and Toczylowski 3 using all-electron double-excitation 
coupled-cluster theory with a non-iterative perturbational treatment of triple excitations (CCSD(T)) and an av5z+ 
Gaussian basis set; and (hi) a fit of the potential of Korona et al? 1 to our DMC energy data, as described in Sec. IIVI 
We believe the HFD-B pair potential to be the most accurate neon pair potential in the literature to date. 



2. Static-lattice energy-volume curve using pair potentials 

Let the pair potential between two neon atoms at R and R' be </>(|R — R-'D- Let A be a large radius. We evaluate 
the static-lattice energy per atom as 

E Sh *l( £ m\) + ^r r^r 2 cf>(r)dr), (4) 

1 Vo<|R|<A V JA ) 

where the {R} are the lattice sites and N-^ e /V is the number density of neon atoms. This expression becomes exact 
as A goes to infinity. The integral in Eq. (j4]) was evaluated analytically for each pair potential, while the sum was 
evaluated by brute force. A was increased until E SL converged. 



3. Force- constant calculations using pair potentials 

Pair potentials were used to generate force-constant data for a quasiharmonio^ calculation of the ZPE of neon. 
The method of finite displacements^ was used to generate the force constants in a finite supercell subject to periodic 
boundary conditions. It was ensured that the force constants were highly converged with respect to the size of the 
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FIG. 2: (Color online) Same as Fig. [TJ but at higher densities. 



displacements and the number of periodic images of the neon atoms that contributed to the force constants. Following 
the evaluation of the force constants, the calculation of the ZPE proceeded as described in Sec. Ill A 21 



D. EOS models 



Let E(V) be the total energy of a neon crystal as a function of primitive-cell volume V. It has previously been 
noticed 32 that a Vinet EOS of the form 



1/3N 



exp 2^o-l) 1 




(5) 



where the zero-pressure volume Vo, bulk modulus Bo, pressure-derivative of the bulk modulus B' , and integration 
constant C are fitting parameters, gives a better fit than a third-order Birch-Murnaghan EOS, 
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V r 



7/3 



E(V) = --B I (4 - B' )^ (14 - SBJ)^ 



(16 - 3B' 



5/3' 



y2/3 



c, 



(6) 



to DFT results for solid neon. In some cases the Vinet EOS gives a lower \ 2 value when fitted to our DMC data than 
the Birch-Murnaghan EOS; in others it gives a higher \ 2 value. For example, using DMC data obtained in simulation 
cells consisting of 2 x 2 x 2 primitive cells and a time step of 0.01 a.u., the Vinet and Birch-Murnaghan EOS models 
give x 2 values of 10.0914 and 35.3561, respectively, whereas at a time step of 0.0025 a.u. the EOS models give x 2 
values of 20.6609 and 4.8967, respectively. The resulting pressure-volume curves are essentially indistinguishable in 
each case, however. To be consistent, we have fitted Vinet EOS's to all of our theoretical data. 



III. DFT STUDY OF PHASE STABILITY AND BAND GAP 



A. Phase transitions in solid neon 



We have compared the DFT energies of face-centered cubic (FCC) and hexagonal close-packed (HCP) phases of 
solid neon. For HCP neon the lattice-parameter ratio c/a was optimized, but the optimal ratio always turned out to 
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Primitive-cell volume (a.u.) 

FIG. 3: (Color online) Low-density static-lattice pressure obtained by fitting Vinet EOS's to the DMC energies obtained in 
simulation supercells consisting of n x n x n primitive unit cells and at different time steps r. The pressure extrapolated to 
infinite system size is also shown. 
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FIG. 4: (Color online) Same as Fig. (3j but at higher densities. 



be y/8/3, which is the ratio appropriate for an ideal HCP lattice. The DFT energy difference between the FCC and 
HCP phases is typically less than 0.0005 a.u.: too small for us reliably to identify any phase transition. Experimentally, 
Hcmley et alM have found that solid neon adopts the FCC phase up to pressures of at least 110 GPa at 300 K. We 
have therefore used the FCC lattice in all of our calculations, apart from those described in this section. 
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FIG. 5: DMC energy of solid neon plotted against time step for a simulation supercell consisting of 3 x 3 x 3 FCC primitive 
unit cells. The primitive-cell volume is 85.75 a.u. 960 configurations were used in the DMC simulations. 



B. Band gap of solid neon 

The band gap of solid neon, calculated using DFT, is shown in Fig. [6l The band gap is large at the equilibrium 
volume, and increases significantly when the material is compressed. The DFT calculations predict that neon is still 
an insulator when it is compressed to a primitive-cell volume of 2 a.u., corresponding to a pressure of about 366 TPa. 
The use of the ultrasoft neon pseudopotential (with a core radius of 0.9 a.u.) probably causes the DFT results to 
become unreliable at such high densities; nevertheless, our results indicate that the metalization pressure of neon is 
of the order of hundreds of TPa. Hemley et aL— concluded that neon remains a wide-gap insulator over the range of 
pressures that they studied using diamond-anvil cells (up to 110.4 GPa), while Hawke et alM- used a magnetic-flux 
compression device to show that solid neon remains an insulator up to at least 500 GPa. 

The DFT-LDA and DFT-PBE band gaps at the experimental equilibrium primitive-cell volume (150 a.u.) are 
11.85 eV and 12.04 eV, respectively, compared with the experimentally determined value of 21.51 eVi^ As usual, DFT 
substantially underestimates the band gap. The GW band gap of 20.04 eV, calculated by Galamic-Mulaomerovic 
and Patterson, is relatively accurate^ The DMC method can also be used to perform highly accurate band-gap 
calculations , 37 ' 38 although we have not done this for neon. 



IV. DMC-CALCULATED PAIR POTENTIAL FOR NEON 



The difference between the DMC pair potential, evaluated as the fixed-nucleus total energy of a neon dimer, and 
the HFD-B potential is shown in Fig. [7] The DMC energy data have been offset by a constant that was determined 
by fitting the data to a pair-potential model.— We have used the form of potential proposed by Korona et al. for 
helhmij^i which can be written as 

4>{r) = Aexp (~ar + pr 2 ) + £ f 2n (r, b)^, (7) 

n— 3 

where r is the separation of the neon atoms. The dispersion coefficients C2n are taken from Cybulski and Toczylowski 3 
(C 6 = 6.28174, C 8 = 90.0503, C w = 1679.45, C 12 = 4.18967 x 10 4 , C u = 1.36298 x 10 6 , and C 16 = 5.62906 x 10 7 ) 
and f2n( r i b) is the damping function proposed by Tang and Toennies^ 

/ 2n (r,6) = l-exp(-6r)X;^. (8) 

k=0 
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FIG. 6: DFT band gap of FCC neon against the primitive-cell volume, calculated using the LDA and PBE exchange-correlation 
functionals. 
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FIG. 7: (Color online) Neon pair potential, calculated using DMC. The statistical error bars on the DMC results are smaller 
than the symbols. A fit of the form of pair potential proposed by Korona et a/.— to the DMC data is also shown, as is the pair 
potential generated by Cybulski and Toczylowski using CCSD(T) theory.— All are plotted relative to the HFD-B pair potential 
of Aziz and Slaman.— 



A, a, f3, and b are adjustable parameters, which were determined by a x 2 fit to the DMC data, as was the constant 
offset. The fitted parameter values (in a.u.) are A = 84.956788, a = 2.0683266, /3 = -0.II767673, and b = 2.6899868. 
The difference of the resulting pair potential with the HFD-B potential is shown in Fig. \7\ as is the corresponding 
curve for the CCSD(T) data. It can be seen that over a wide range of separations the DMC pair potential lies closer 
to the HFD-B pair potential than the CCSD(T)-generated pair potential. 
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FIG. 8: (Color online) Phonon dispersion curves calculated using DFT and pair potentials for FCC neon at a primitive- 
cell volume of 41.59375 a.u. The Einstein frequencies evaluated using DFT-LDA, DFT-PBE, the HFD-B pair potential, the 
CCSD(T) pair potential, and the DMC pair potential are 83.10234 meV, 87.57110 meV, 63.77843 meV, 65.73911 meV, and 
63.00858 meV, respectively. The quasiharmonic ZPE evaluated using the HFD-B potential is 0.003321155 a.u. 



V. LATTICE DYNAMICS AND ZERO-POINT ENERGY 

By comparing the results obtained using a Lennard-Jones potential in the harmonic approximation with the VMC— 
results obtained using the same potential by Hansen,— Pollock et aL— have demonstrated that the harmonic approx- 
imation is valid for solid neon at high pressures. We consider two methods for calculating the ZPE: (i) the ZPE of 
quasiharmonic phonons can be evaluated in a supercell of several primitive cells, or (ii) the ZPE can be computed 
within the Einstein approximation by evaluating the quadratic potential felt by each atom as it is displaced from its 
equilibrium position with all the other atoms held fixed. 

Examples of phonon dispersion curves at two different densities are shown in Figs. [5] and [51 Inelastic neutron- 
scattering data^ are also shown in Fig. [9] At high density the DFT-LDA and DFT-PBE dispersion curves are in 
good agreement, but at low density the DFT-LDA phonon frequencies are significantly lower than the DFT-PBE 
frequencies. Unstable (imaginary) phonon modes start to occur at a primitive-cell volume of about 133 a.u. in the 
LDA. By contrast, there are no unstable phonon modes, even at a primitive-cell volume of 182.25 a.u., when the 
PBE functional is used. The DFT and pair-potential results are in agreement at high densities, indicating that the 
DFT results are accurate in this regime. Overall, the DFT-PBE dispersion curves appear to be more accurate (that 
is, closer to the HFD-B and experimental results) than the DFT-LDA dispersion curves. 

The pressure arising from the ZPE of solid neon as calculated using different methods is plotted relative to the 
HFD-B results in Fig. [TO] Within DFT, the Einstein approximation is excellent. It can be seen that the difference 
between the LDA and PBE results is appreciable at low densities, but that the difference between the Einstein and 
quasiharmonic zero-point pressures is more significant at high densities. As expected from examination of the disper- 
sion curves, the HFD-B zero-point pressure is closer to the DFT-PBE results than the DFT-LDA ones; nevertheless, 
all the zero-point- pressure results are in good agreement. 

The DFT quasiharmonic zero-point pressures have been added to the corresponding static-lattice pressures to give 
the final EOS's. The DFT-PBE quasiharmonic zero-point pressure has been added to the DMC static-lattice pressure 
to give the final DMC EOS. For the pair potentials, the quasiharmonic zero-point pressure calculated using each pair 
potential has been added to the corresponding static-lattice pressure in order to obtain the final EOS. 
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FIG. 9: (Color online) The same as Fig. [5] but with a primitive-cell volume of 149.06894 a.u. (close to the experimental 
equilibrium density). Experimental data from Ref. are also shown. The imaginary frequencies of unstable modes are 
plotted as negative numbers. The Einstein frequencies evaluated using DFT-PBE, the HFD-B pair potential, the CCSD(T) 
pair potential, and the DMC pair potential are 5.37853 meV, 4.07895 meV, 6.73975 meV, and 6.06113 meV, respectively. The 
quasiharmonic ZPE evaluated using the HFD-B potential is 0.000216603 a.u. 
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FIG. 10: (Color online) Difference of zero-point pressure of FCC neon calculated using various methods and the HFD-B result 
for the zero-point pressure. (The noise is due to the fact that Monte Carlo methods were used to sample the first Brillouin 
zone when calculating the zero-point energy, and the resulting curve was differentiated numerically to obtain the zero-point 
pressure.) 
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FIG. 11: (Color online) EOS of FCC neon, obtained by experiment and various theoretical techniques. 



VI. ZERO-TEMPERATURE EOS OF NEON 



Zero-temperature EOS's for solid neon, calculated using DFT-LDA, DFT-PBE, DMC [extrapolated to infinite 
system size using Eq. and pair potentials are shown in Fig. [TTJ and the differences of the theoretical EOS's with 
the experimental EOS are plotted in Fig. 1121 (The low-density experimental pressure- volume data of Anderson et o/.— 
shown in Fig. [TT] were obtained at 4.2 K, while the high-density experimental data of Hemley et al^ were obtained 
at 300 K. Hemley et al. reduced their pressure-volume data to the zero-temperature isotherm using a Mie-Gruneisen 
model and fitted their results to a third-order Birch-Murnaghan EOS, with the zero-pressure primitive-cell volume 
and bulk modulus being the values obtained by Anderson et al. The resulting EOS 46 is valid at both low and high 
densities and is regarded as being the definitive experimental EOS.) The parameter values for a Vinet fit (Eq. ^) to 
our DMC data (including the DFT-PBE ZPE) are V = 128.52597 a.u., B = 2.7539319 GPa, and B' = 7.6510744. 

At low densities the DFT-LDA and DFT-PBE EOS's differ markedly. The strong dependence of the DFT results 
on the choice of exchange-correlation functional implies that the description of van der Waals bonding within DFT is 
unreliable, as one would expect, given the local nature of the approximations to the exchange-correlation functional. 
DMC produces a considerably more accurate EOS than DFT, suggesting that DMC is capable of giving a proper 
description of van der Waals bonding. The HFD-B pair potential 29 gives an EOS of similar accuracy to the DMC 
EOS at low to intermediate densities. At higher densities, the EOS calculated using DMC is better than any of the 
pair-potential EOS's. Although the difference between the DMC and experimental pressures is significant at high 
densities, it should be emphasized that the fractional error remains small. 

The pair potential calculated using CCSD(T) theory^ gives a poorer EOS than the DMC-generated pair potential. 
On the other hand, the EOS obtained using the DMC pair potential is significantly poorer than the HFD-B EOS. 
Taken together with the fact that the direct DMC EOS is excellent, this suggests that many-body interactions play 
a significant role in solid neon, and that such interactions are included to some extent in the HFD-B potential. 



We have performed DMC calculations of the energy of FCC solid neon as a function of the lattice constant and 
the energy of the neon dimer as a function of atomic separation. Other calculations using DFT methods and pair 
potentials have been performed to evaluate the ZPE and for comparison purposes. 

We have calculated the phonon dispersion curves of solid neon using the DFT-LDA and DFT-PBE methods, the 
HFD-B pair potential, and CCSD(T)- and DMC-derived pair potentials. We believe the results obtained with the 
HFD-B pair potential are likely to be the most accurate. DFT-PBE gives more accurate dispersion curves than DFT- 
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FIG. 12: (Color online) Comparison of theoretical EOS's of FCC neon. The difference of the calculated pressure is plotted 
against the experimental pressure as a function of volume. 



LDA, for which the phonon frequencies are too low. The dispersion curves obtained with the DMC pair potential are 
more accurate than those obtained using either DFT or the CCSD(T) pair potential. We have calculated the ZPE 
of solid neon using the DFT-LDA and DFT-PBE methods, and the HFD-B, CCSD(T), and DMC pair potentials, 
within the quasiharmonic approximation. At low pressures the ZPE depends on the calculation method used, but 
the contribution to the EOS is small, while at high pressures the dependence on the calculation method is relatively 
weak, although the contribution of the ZPE to the EOS is significant. The Einstein model gives ZPE's in very good 
agreement with the quasiharmonic values over the pressure range considered. 

We have calculated the zero-temperature EOS of solid neon using the DFT and DMC methods, including corrections 
for the ZPE. We have shown that the DFT results depend strongly on the choice of exchange-correlation functional, 
while the DMC results are close to the experimental EOS. We therefore have evidence that DMC gives a better 
description of van der Waals bonding in real materials than DFT. At high pressures the DMC EOS is closer to the 
experimental results than the EOS obtained using the HFD-B pair potential. However, the statistical errors of about 
0.0002 a.u. in the DMC energy data for solid neon are too large to determine an accurate value for the lattice constant 
of the solid. We have shown that the neon pair potential determined by DMC calculations gives a more accurate EOS 
than the pair potential determined by CCSD(T) calculations, although the DMC pair-potential results are not as 
accurate as those obtained using the semiempirical HFD-B potential. Overall, our results demonstrate the accuracy 
and reliability of the DMC method and the high quality of the neon pseudopotentials that we have used. 
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